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Abstract 
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A multiphase mathematical model for suspension 
polymerization of VCM in a batch reactor has been presented here. 
This model incorporates mass transfer for the initiator, monomer 
and radical species between the different phases involved Ci.e. , 
monomer, polymer, VCM droplet surface, aqueous and vapor phasesD. 
This includes the possibility of inhomogeneity of both the 
initiator distribution and initiator efficiency. The model 
accounts for precipitation of polymer from monomer phase and also 
the pressure drop in the reactor. The model has been modified to 
incorporate the continuous addition of seal and Injection water 
in accordance with industrial practice. The model predicts the 
rate of polymerization, pressure drop in the reactor, 
instantaneous and cumulative molecular weight averages, short and 
long chain branching characteristics with respect to process 
conditions. The sensitivity of polymerization rate and pressure 
to various model parameters has been investigated. These 
sensitive parameters have been tuned with respect to industrial 
data. The parameters have been estimated using the simplex 
optimization technique. When the objective is to utilise the 
maximum cooling capacity of reactor the optimal temperature 
profile has been obtained by using Pontryagin’s Prlnci|>le. We 
have also estimated various heat transfer coefficients using 


industrial data. 
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CHAPTER 1 
INTRODUCTION 

1.1 Overview 

Polyvinyl chloride CPVCD is one of the oldest thermoplastic 
polymers. The commercial production of PVC started in 1631 in 
Germany and it became commercially very significant prior to 
World War II. A major reason for the long delay between discovery 
and commercialisation was that unmodified PVC is not a useful 
material. A significant step in promoting the growth of PVC 
industry was the development of o5<ychl or i nation process which 
made it possible to use ethylene as feed stock for vinyl chloride 
monomer CVCMD production. Another major factor for the growth of 
PVC industry is that it is the largest consumer of chlorine, which 
is essentially a by-product in the manufacture of sodium 
hydroxide. Currently, next to polyethylene, PVC is the most 
frequently used thermoplastic material. 

The commercial production of PVC is based on three 
processes: suspension, emulsion and bulk polymerization. The 
suspension process is the dominating route to PVC. It accounted 
for about 80$^ of world production in 1680. In India 
also, suspension process accounts for the major production 
capacity. The industries that manufacture PVC in India at present 
C1663!> are Indian Petrochemicals Corporation Ltd, Chemicals and 
Plastics India Ltd, Reliance, Finolex, Shri ram Plastics, etc.. 

The suspension process is favoured because of the relative 
ease of polymerization control Cdue to water acting as a large 
heat sink3 and lower cost of separating and drying the resin. The 
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suspension process is carried out in a glass lined or highly 
polished st^ainless s'teel st-irred tank batch reactor with a 
jacket, meant for heating and cooling. The raw materials used are 
Vinyl Chloride Monomer CVCMD , demineralised water and certain 
additives including a suspension stabiliser. Polymerization of 
VCM is a free radical type of chain reaction. Once the batch 
attains the desired reaction temperature by heating, the 
polymerization reaction becomes self sustaining. It is exothermic 
in nature and heat has to be removed continuously for stable 
operation. The heat of reaction is removed through the jacket to 
bring the reactor contents to the desired polymerization 
temperature. Normally it takes 3-6 hrs to reach 85-00% 
conversion. The process is characterised by an auto-accelerating 
real time profile ,in some cases producing a marked peak near the 
end of the reaction, which can cause problems with temperature 
control and the danger of run away reaction. 

The product quality of suspension PVC is, to a large extent, 
controlled by the particular suspension stabilizer system used 
and the precise manner in which the polymerization is carried 
out. The suspension stabilizer system affects the shape of the 
grains, the grain size distribution, resin porosity and bulk 
density of resin. Typical suspension stabilizers used in industry 
are methyl hydroxy propyl cellulose, poly vinyl alcohol etc. 

1.2 Literature Survey 

There has been considerable amount of work on the modelling 
Cl-lSl of PVC formation reported in the literature, particularly 
in developing kinetic mechanisms. These models are all based on 
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the key feature that PVC is essentially insoluble in its monomer 
so that as polymerization proceeds within the monomer 
dr opl ets , sol i d polymer precipitates out and swells to the extent 
of approximately 23J^ by weight of VCM to form a deformable 
vi scoel asti c gel . 

1 . 2. 1 Two Phase Model 

Talamini Cll were the first to derive a two phase C monomer 
rich phase and polymer rich phaseD model to predict the time 
variation of fractional monomer conversion in a suspension PVC 
batch reactor. The assumptions made in Talamihi*s model are 
i3 The corresponding initiator decomposition and propagation 

rate constants as well as the initiator efficiencies in the 

two phases are equal 

ii3 The initiator concentration is the same in both the phases. 
iiiZ) No transfer of radicals between two phases occurs. 

Iv5 The initiator concentration does not change with conversion. 
v3 The volume of the reaction mixture remains constant. 

This was the first model to predict the rate of polymerization up 
to 70Ji conversion. The onset of two separate phases was reported 
to begin after less than conversion and lasts until between 
70-805^ conversion. This particular conversion at which the 
monomer phase gets depleted is called the critical conversion x^. 

Abdel -Alim and Hand elec [21 modified Talamini's model by 
replacing assumptions Civ5 and CvD with the following 
assumptions: 

viD The initiator concentration changes with reaction time. 
viiD The volume of the reaction mixture varies linearly with 
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conversion. 

This model assumes the presence of two phases in equilibrium»each 
of constant composition Cx^& of polymer respecti velyD . x^ is 

very small, so taken as zero, that is, the monomer rich phase 
contains only monomer. This model could also predict the 
molecular weight distribution CMWDD throughout the entire range 
of conversions and has been applied to commercial PVC processes. 

Ugelstad 131 assumed production of radicals in both phases 
in contrast to Abdel -Alim and Hand elec model. They also included 
the transfer of radicals between the two phases. However, they 
suggested that an equilibrium distribution of radicals is quickly 
established between the two phases, which means no net transfer of 
radicals takes place between the phases. Ugelstad [41 and OlaJ 
151 have put considerable efforts into refining the kinetic 
mechanism of the process, but have tended to concentrate on the 
low conversion behaviour. Kuchanov and Bort 161 assumed that the 
radical transfer from polymer phase can be completely neglected. 
Whilst these earlier models recognised the existence of the two 
distinct phases in the polymerization, they do not consider the 
detailed physical nature and structure of these phases. 

1.2.2 Multiphase Model 

Kelsall and Maitland 171 proposed a multiphase mass transfer 
model Cmonomer, polymer, VCM droplet surface, aqueous, and vapor 
phases!) , incorporating: CaZ) mass transfer for various species 
within the different phases involved, and CbD inhomogeneous 
initiator distribution between monomer and polymer phase and CcZ) 
Initiator efficiency ratio between monomer and polymer phase. The 
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model can predict rate profile, and several important properties 
of product like porosity, apparent density and particle size etc. 
from a knowledge of process conditions. But this model does not 
account for diffusion controlled propagation , termination rate 
constants at high monomer conversions and transfer to polymer. 
1*2.3 At High Conversions Cbeyond 

Many workers [8-151 have studied experimentally the 

polymerization of vinyl chloride. Xie [81 experimentally 

investigated the vinyl chloride polymerization at high conversion 
and reported that reactor pressure falls before monomer is 

completely consumed and the conversion at which pressure drops 
depends on the morphology of the PVC particles. They have 
surmised that reactor pressure falls in two stages. The first 
stage is due to the volume increase of the vapor phase as a 
result of volume shrinkage during reaction. The monomer phase has 
not yet been consumed at this stage but is trapped in the 

interstices between primary particles creating a mass transfer 
resistance; therefore the pressure drops slowly. The second stage 
is due to both the volume increase of vapor phase and to the 
monomer in the vapor phase diffusing into polymer phase. They 
have said that the thermal stability and porosity of PVC. product 
decreases significantly with conversion beyond the critical 
conversion. 

Xie [91 proposed an equilibrium model relating 
temperature, pressure, monomer conversion, and monomer phase 
distribution for vinyl chloride polymerization. This model can be 
used to determine the monomer conversion beyond the pressure drop 
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by measurement of reactor temperature and pressure. Various 
assumptions have been incorporated into this model like vapor 
phase inside the reactor follows ideal gas law and the solubility 
of VCM in water follows Henry’s law. 

Hand elec [lOJ first applied the free volume theory to model 
diffusion-controlled termination and propagation to explain the 
properties of PVC obtained at high conversions, which has been 
utilised by Xie 1111. These modelling attempts were premature due 
to inadequate amount of kinetic data at these high conversion 
levels. 

1.2.4 Semi -batch Operation 

Xie 1121 has developed a semi -batch reactor model for 
suspension polymerization of vinyl chloride. In the conventional 
batch process, at high conversions Cconversion > D, the 

kinetic behaviour of VCM polymerization changes with conversion 
dramatically, and this includes changes in reactor pressure, 
monomer concentration, reactor rate, instantaneous and 
accumulated molecular weight averages. Also, properties of PVC, 
such as thermal stability, deteriorate with increasing 
conversion. All the kinetic features like decrease in 
decomposition rate constants, propagation rate constants, 
efficiency factor, and increase in chain transfer to monomer rate 
constants which occurs at high conversions can be attributed to 
the decrease in monomer concentration. He has found 
experimentally that these changes are minimized if a constant 
monomer concentration is maintained at high conversions by 
semi— batch operation. Also thermal instability of PVC is found to 
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be minimized at high conversions by keeping monomer concenlrabion 
constant.. 

1.2.5 Polymer Properties and Thermal Stability 

The effect of polymerization conditions on polymer 
properties have been investigated by the same authors 1131 and 
they have found that thermal stability of PVC decreases with 
increase in polymerization temperature and the reactions forming 
defect structures are favoured at low monomer concentrations. 
Thermal instability of PVC is attributed to the formation of 
allylic and tertiary chlorine which becomes prominent after 
x^.They have experimentally found out that dehydrochlorination 
rate increases as conversion increases beyond x^» but it is 
independent of conversion below x^. It also increases as 
temperature is increased. They have suggested, therefore, the 
most effective way to minimize the formation of defective 
structures is to maintain monomer concentration as high as 
possible during polymerization after pressure drop. 

The most recent study on modelling of PVC suspension process 
was done by Sidiropolou and Kiparissides 1141, where they used a 
similar approach as used by Kelsall and Maitland [71. This model 
can give some new results like, the number of short and long chain 
branching as well as the number of unsaturated terminal double 
bonds per polymer molecule. This new approach may be useful for 
commercial PVC production. 

Xie (111 has studied the kinetics and mechanism of vinyl 
chloride polymerization in detail from both chemical and physical 
point of view. He has said that chain transfer to monomer plays a 
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vital role in controlling molecular weight development in VCM 
polymerization. Also, the limiting size of the primary particle 
is approximately 1.4pm in diameter and the final number of 
primary particles is around 2.0*10^^ cm~^ but is a function of 
reaction temperatures. 

The molecular weight development during the suspension 
polymerization of VCM has been investigated by Xie [151. He 
states that at high conversions the number average molecular 
weight decreases significantly as compared to weight average 
molecular weight consequently increasing the polydispersity 
index. The molecular weight of PVC increases significantly with 
decrease in polymerization temperature. 

1 . 2. 6 Optimization 

Albright and Soni 116] have investigated the suspension PVC 
p>olymerlzation with multiple initiators. They say by using a 
mixture of fast and slow initiators, it is possible to get a 
completely flat rate profile. Economic saving of at least 10 - 
30^ can be realized by changing to initiators causing flatter 
rate of heat release. Large economic savings are also possible if 
a reflux condenser can be used to remove a portion of heat of 
reaction and of agitation from the reactor. The same authors I17J 
have developed a computer program for determining the optimum 
design of suspension PVC batch reaction system. Their model 
calculates design and key operating variables, number of reactors 
required for given annual production rates, capital costs, and 
operating costs. 
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1.2.7 Diffuslvlty 

Berens 1181 has i nvest-i ga'ted "the diffusion process of VCM in 
PVC. He has found by experimen'ts 'Lha't initially for short time 
fickian diffusion takes place and later on the relaxation by 
swelling dominates thus producing non - fickian diffusion. 

1.3 Mor phol ogy of PVC 

The suspension process produces a particulate, free flowing 
product having the same general appearance with a mean particle 
size usually in the range 100-150/jm. 

Suspension polymers are usually porous in nature and made up 
of a large number of much smaller particles. Suspension PVC 
particles usually possess a pericellular skin or membrane which 
extends almost continuously over the entire outer surface of 
particle. This skin is known 1101 to form during the early stages 
of polymerization C < Sti conversion D and this has been shown to 
be a PVAyPVC graft copolymer . It has been postulated 1201 that 
the skin is formed by primary particles that, when first 
generated, are brought to monomer /water interface of dispersed 
droplets by centrifugal forces. They then destabilise, such that 
they fuse and form a continuous boundary that eventually grows to 
a thickness of about l/jtn. As conversion rises the skin- becomes 
rigid depending on the recipe used. Varying levels of coalescence 
within the skin now take place over the range of conversion 6 - 
1B>4 1211. At a conversion close to critical conversion, the PVC 
Internal structure and the skin surrounding the particles becomes 
rigid enough to give a significant resistance for monomer 
diffusion from vapor phase to polymer phase [81. 
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The morphological properties of the suspension PVC grains 
include shape, size & its distribution, surface characteristics 
internal structures and porosity. Among them, the size and its 
distribution, as well as porosity have crucial influence on the 
separation of free unreacted VCM from the grains and on the 
processing and the end use of these resins. 

Process Variables Affecting the Particle Properties of Suspension 
PVC 

1 . 3. 1 Polymerization Tenperature 

In the absence of other reactive agents, the molecular 
weight of PVC is almost entirely determined by polymerization 
temperature 1191. Increasing the polymerization temperature in a 
reactor reduces the colloidal stability of domains and causes 
earlier aggregation. The size of primary particles increases and 
the number decreases, which leads to the formation of weaker 
networks of aggregates, less able to resist droplet contraction, 
and consequently grains with lower porosity are produced. 

1.3.2 Conversion 

Conversion is found to have a very large effect on product 
properties. In the range of conversions commonly encountered 
commercially Ci.e. , 75— OSJi > porosity is generally found to 

decrease with increasing conversion. This is because very low 
conversion products, lacking in internal strength collapse 
inwardly on themselves thus reducing porosity and increasing bulk 
density. 
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1.3.3 Dispersant 

The polymerising monomer droplets in the suspension process 
are formed and stabilised by the combined effect of vigorous 
agitation and the presence of a dispersant or dispersants. The 
choice of the dispersant is of utmost importance as it not only 
contributes to controlling the particle size distribution of PVC 
produced but also has a major effect on the sub-structure of the 
particle in terms of its porosity etc. 

1.3.4 Agitation 

Agitation is of fundamental importance in the suspension PVC 
process. Together with the dispersant system it governs the 
stability of the suspension during polymerization and the 
particle size of the product formed. Agitation can also have an 
important role in determining other product properties such as 
porosity and bulk density. 

In an agitated dispersion, the turbulent velocity 
fluctuations through out the mixing tank tend to bring about 
separation of all adhering droplet clusters during the critical 
time period, so that the chance of coalescence is minimal. Both 
the force of adhesion between droplets and the forces caused by 
agitation C which tend to separate the dropletsD depend on the 
diameter of the droplet. However, the changes of separation are 
greater the larger the droplets and thus, for a given level of 
agitation, there exists a minimum droplet diameter above which 
stabilisation by agitation becomes possible. By contrast, the 
susceptibility towards droplet break-up, caused by local velocity 
fluctuations and shear forces near the impeller, is increased by 
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increases in droplet size. The minimum droplet size that is able 
to remain stable to break-up will, of course, increase with 
decreasing agitator speed. However, in the case of very low 
agitation intensity large droplets may become unstable due to 
phase separation cause by specific gravity differences. 

Increased severity of agitation is often found to increase 
porosity although this does not appear to always be the case. 
Hofmann [221 in fact found an initial reduction in porosity with 
Increasing agitator speed followed by an increase corresponding 
to the increase in particle size observed. He attributes the 
increase in product porosity that occurred under high agitation 
conditions to two factors. Firstly, the increased degree of 
agglomeration results in the formation of a greater number of 
voids at the interstices between the discrete particles. 
Secondly, increased porosity is said to occur due to the internal 
structure of the primary particle being disturbed, with the 
membranes of the individual polymerising droplets being torn and 
partially destroyed by the agglomeration process. 

1.3.5 Other variables 

Water : monomer phase ratio is often considered an important 
production parameter in suspension polymerization. Certainly in 
the case of vinyl chloride polymerization there are physical 
limitations to the minimum phase ratio. A very low proportion of 
water will cause instability and gross distortion of the particle 
size distribution. At a level somewhat higher than this, 
depending on resin porosity, problems are encountered due to 
absorption of water within the grains causing a drying out effect 
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and thus high slurry viscosities. This effect can be overcome by 
adding extra water during polymerization. 

The presence of even trace quantities of oxygen during the 
polymerization can have dramatic effects on both the chemical and 
physical properties of the resulting PVC. The most obvious effect 
is that on particle size which is decreased considerably as the 
reactor oxygen content is increased. 

Cebollada [231 has found by experiments that molecular 
weight of PVC remains insensitive to changes in viscosity of the 
suspension medium. For all practical purposes* viscosity bears no 
influence on the reaction kinetics. But it plays a fundamental 
role in the determination of particle size distribution. High 
viscosity media produce larger size* unicellular spherical 
particles retaining their identity as individual droplets. 
Conversely* low viscosity media favor the formation of smaller 
particles. 

The processing character sties of suspension PVC are largely 
determined by the molecular weight of the polymer and morphology 
of the grains (241. As molecular weight is reduced* the viscosity 
of the hot melt decreases and polymer becomes easier to process, 
although the strength of the finished article is reduced. 

1.4 Proposed Objectives 

The objectives of the present work are the following; 
aD To simulate the multiphase mass transfer model for 
suspension PVC polymerization process in a batch reactor. 
b3 To find out what model parameters are sensitive to rate of 
polymerization and pressure drop in the reactor. 
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cD To estimate© "the aforement-ioned parameters using 

industrial data. 

dD To estimate the heat transfer coefficients from the 
industrial data. 

eD To find out the optimal temperature policy for 

maximum utilisation of cooling capacity of the reactor. 

fD To predict instantaneous and cumulative molecular weight 
aver ages . 

g5 To use semi— batch mode of operation with respect to monomer 
after critical conversion to keep the pressure constant, 
reduce the thermal instability, and to increase the 


productivity. 
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CHAPTER 2 

MULTI -PHASE MODEL DEVELOPMENT 
2. 1 The Suspension Process 

Before giving ■the derails of t-he model, let us see the 
salient features of suspension PVC process it attempts to 
simulate. 

3 

The reactor considered here is a 70m reactor. The initiator 
used is 2-ethyl hexyl per oxydi carbonate CEHPD. VCM and DM water 
Cdemi neral i sed waterD are added simultaneously into the reactor. 
EHP is then added over a short time. During polymerization seal 
water is added continuously in to the reactor to prevent the back 
flow and injection water is added to maintain the liquid volume 
constant. The excess heat from the reactor is removed by two 
ways. Coolant water is circulated through the half pipe jacket 
and the chilled water is circulated through the baffles inside 
the reactor. There are four baffles in the reactor with two 
concentric pipes where coolant flows from bottom to top in the 
inner pip© and drains out in the outer pipe. The schematic 
diagram of the reactors are shown in Fig 2.1a and Fig 2.1b. 

PVC polymerization can be visualised in two ways: one is 
macroscopic scale and the other is the microscopic scale.' 

Macroscopic Scale : An aqueous suspension of VCM droplets at 
a volume fraction of approximately 40 is formed in the autoclave 
by vigorous agitation. The initial droplets are typically 
50-150/./m in diameter and are stabilised by the addition of a 
protective coalescence. Depending on the nature of the agitation 
and the degree of stabilisation, the pol 3 anerizing droplets may 
aggregate to a certain extent as conversion increases to produce 



12 



Figr 2. la Schemat^ic Represent.a'tion of" Ihe Reactor 


1 VCM charge 

2 Jacket coolant inlet 

3 Baffle coolant outlet 

4 Baf f 1 e cool ant 1 nl et 

5 Seal water addition 

6 DM water charge 


7 Agitator 

8 Injection water addition 

9 Baffles 

10 Half pipe Jacket coil 

11 Jacket coolant outlet 

12 Initiator charge 






a final particle C grain!) of approximately 50-1 SOpm. In 
addition, agitation must give adequate heat transfer from the 
droplets to ensure uniform temperature control. 

Microscopic Scale : Superimposed upon this droplet 

aggregation pr ocess , another sequence takes place within each 
droplet. The growing chain radicals become insoluble in VCM at 
relatively short chain lengths C appr oxi matel y 10 monomer unitsD 
and continue to grow as precipitated coiled macro-radicals. The 
first chains so formed aggregate in units of about 50 to produce 
small particles called micro-domains of diameter 200A°. These in 
turn flocculate to produce domains or primary particle nuclei 
CO. 1 to O. 2piiD by 1?< conversion. All growing macro-radicals and 
dead piolymer chains are now absorbed by these units before they 
have a chance to nucleate fresh micro-domains and for the period 
the total number of domain particles remains constant, and they 
simply grow in size. These primary particles CO. 2-0. 4pm!) 
eventually become unstable and flocculate, from about 3?< 
conversion onwards, to form close-packed agglomerates of l-2pm 
diameter. The way in which these units grow and pack together 
within the droplet is crucial for the development of both the 
reaction rate,grain porosity and apparent density. These primary 
particle aggregates continue to grow up to 10pm in diameter by 
surface deposition and internal growth, with increasing 
intra-primary growth, until polymerization is terminated at 
typically .90-955^ conversion. This sequence of microscopic growth, 
which determines the essential char acterstics of the final 
product, is summarised schematically in Fig 2.2. 
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rig 2.2 Scheme of PVC 
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2. 2 Nat/Ure of Phases 

A realist-ic model of t,he process should include 

polymerizat-ion wit,hin and mass t-ransfer bet, ween all t,he phases. 
There are five phases » viz, monomer rich phase C phase ID, polymer 
rich phase C phase 2D, VCM droplet, surface phase C phase 3D, 
aqueous phase C phase 4D , and vapor phase C phase 5D. 

The major loci of polymerization are the two immiscible 
PVC/VCM phases. 

Monomer Rich Phase 

This is essentially pure monomer at conversions, less than 
0. 001;this is of course, the sole site of polymerization. This 
phase disappears at conversions approximately equal to 

0.77, beyond which there is insufficient monomer to fully swell 
the formed polymer in polymer phase, and the autoclave pressure 
falls. The initiator is chosen to be highly soluble in this 
phase; chain radicals initiated here will precipitate out as soon 
as they grow more than a few monomer units into polymer phase. 
These precipitated chains are designated as phase 1 chains as 
long as they continue to grow, transfer to monomer or terminate by 
combination/disproportionation with another radical whilst 

suspended in phase 1. On adsorption by the phase 2 primary 
particles they become phase 2 radicals. 

Polymer Rich Phasd 

This has an equilibrium composition of 0.77 mass fraction of 
polymer and 0.23 mass fraction of monomer and exists as highly 
deformable primary particles which grow and coalesce in manner 
described in section 2.1. The high density of chains in this 
phase has the effect of significantly reducing the rjate at which 
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the end segments of the growing chain radicals can diffuse 
towards one another and couple to produce a dead chain. The 
simplest way to model this is an overall reduction in the 
effective rate constant for chain termination. The initiator is 
usually highly soluble in this phase also, but not necessarily to 
the same extent as in the phase 1. The creation of dead polymer 
is therefore accompanied by the transfer of both monomer and 
initiator from phase 1. Polymerization in phase 2 therefore 
occurs through Ca5 chains which were initiated in phase 1 and 
have migrated in to the phase 2 prior to termination or transfer 
and CbD chains which have been initiated in phase 2. There is 
evidence I25J that the rate of initiation process can be 
significantly lower than that in the phase 1. 

VCM Droplet Surface 

Electron microscopy studies 1261 indicate that the surface 
region or skin of the initial subgrains is significantly 
different from the bulk in that the number density of primary 
particles is higher here, leading to an enhanced rate of 
polymerization in the early stage of process. In addition, the 
chemical nature of the polymer formed in this surface region can 
vary by incorporation of the adsorbed surfactant to form a graft 
copolymer. Both these factors contribute to the subgrain surface 
attaining a coherent close packed structure far earlier in the 
polymerization than the bulk. 

Aqueous Phase 

This phase fulfils the major role of providing a large heat 
sink through which the heat released by pol ymer i zati on may be 
transfered to the autoclave wall. Although the sol ubi 1 i ti es df 
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VCM and initiator in water are both small, they may be significant 
in a full treatment , both in terms of polymerization and 
parti cularl y, in affecting the rates of mass transfer of these 
species in to the polymerizing grains. 

Vapor Phase 

A few percent of the VCM in the autoclave resides in the 
vapor phase. This mass increases slightly as polymerization 
proceeds due to the volume decrease of the condensed phase. After 
pressure drop this remains as the sole reservoir of free monomer 
available to swell the gel phase; the rate of transfer of this 
monomer through phase 4 to the polymerization zone can greatly 
influence the rate of pressure drop. 

2. 3 Kinetic Mechanism 

The kinetic mechanism of VCM polymerization has been 
investigated in much detail, and comprises many side reactions. 
In this work a simplified mechanism has been used, and the 
kinetic scheme is given in Table 2.1. The various options of 
reactions and mass transfers taking place in different phases is 

shown in Fig 2.3. 

/ 

2* 3. 1 Kinetic Parameters 

The kinetic parameters used in our model is given in 
Table 2.2. Although Hamielec tlO] has reported that initiator 
decomposition rate constants, initiator efficiency, and chain 

transfer to monomer decreases with conversion after x , when we 

c 

used his rate constant equations our model did not produce 
satisfactory results. It produced a steep fall in polymerization 
rate at high conversions. 



Fig 2.3 Options for reaction and mass transfer within and between 
different phases 
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C Pol ymer Phas eD 
decomposi ti on , pr opagati on 
termination and chain transfer 
dead polymer present 



CSkin Phased 
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termination and chain transfer 
dead polymer present 

- , ■ 
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Table 2.1 


Kinetic Scheme for Vinyl Chloride Polymerization 


Initiation: 


T ’^di 

I2 *' 


r . k CI_] 

di 1 di 2 


Pr opagati on : 


R + M 
n 


r . = k . [R ] [M3 
pi pi n 


Chain transfer to monomer: 


R + M 
n 


'trMi 


M 

^tri 


k. [R 3 tM3 
trMi n 


Chain transfer to polymer: 


R + P 
n m 


trPi 


P + R 
n m 


k. CR 3 [P 3 
trPi n m 


Ter mi nati on : 


R + R 
n m 


P + P 
n m 

Cor P ^ D 
n+m 


k. . CR 3 [R 3 
ti n m 


i refers to the i phase.i =1 , 2, . . . ,5 


25 


TABLE 2. 2 

Kinetic Rate Constants Used in the Model 


Two options for using various kinetic rate constants are 
below. 


15 Decomposition rate constants 


k .. = 1 . 5»tl0 

dl 

k = k 
d2 '^dl 


15 


expC -14554. 0/15 


at any x 


Cor 5 


*^d2 ” *^dl 

k^g “ *^cjl Cl.O/V 


fp 


- i.O/V 


fx 


55 


X < X 


X > X 


where C = 477.0 expC -2201/10 
ii5 Propagation rate constants 


pi 


5.0»tl0 expC -3320/15 


k o = k , 
p2 pi 

’'p2 = •'pi <-0.2C1.0/V^p 


1.0/V. 55 

fx 


X < X 


X > X 


Cor 5 


k _ = k , exp <-B Cl. O/V^ - 1 . 0/V. 55 

p2 pi ^ fp fx^ 

M 1 

where B = 1.85*tl0 expC -2505/15 


X > X 


iii5 Termination rate constants 
1 P 

k^j^ = 1.3>tl0 expC-2100/D 


t 20 


kt^/a 


2. O 


k = k <C1 
t2 t20 


N yr* 0 - .^2.0 

X 5/Cl -x5> Ck _/k . 5 
c p2 pi 


X < X 


X > X 


where ot 


27.0 - 0.14 CT 


2735 


Cor 5 


’^ta = ’^ti ""'^P - "^cr"> 

where A** = 6. 64*HO® expC -4086/15 


at any x 


given 

Ref 


[14] 


[111 


1141 


rioi 


1141 


[10] 
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i vD Chain tranfer rate constants 
For monomer , 

ktrMi = 5.78 exp C-2768/T:) 

^trMa = ’"trMl ^ 

Cor5 

^trMl ^ 2863.5 exp C-4262. 0/TD 

^trM2 827.11 exp C-4476. 83/10 k^^ x: ^ 


^trM2 




X > 


where K. = CK. D exp <-CH Cl/V. - l/V. D> 

1 1 X fp fx 

c ^ c 

CH** = 1.28»tl0^ exp C -3250/D 

= 5.55*10^^ exp C -9340/D 
c 

CK./KD = CK./KD exp-CC -B*C1 /V. - 1/V. 5> 

4 5 4 5 X ^ fp fx 

c ^ c 

CK-/K-D = 4260 exp C-3410xTD 

4* O X! 

C 

For pol ymer » 


X 

c 


32 Jc 1 i-^ 

■trPl -trHl 

k = k 

^trP2 trPl 

CorD 

kj. r !.4 = 8. 31^10® exp C-11100/T> k . 

t^rPi pi 

^trP2 “ ’^trPl 

v3 Initiator efficiency factor 


at any x 


at any x 


f. =1.0 
1 

f = 0. 1 - 1 . O ' at any x 

£4 

CorD 

f. = 1.0 

1 : ^ 

fg = 0.1 - 1.0 X < x^ 


[141 


f 151 


[141 


1151 


171 
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‘‘a = '-Bf 


1/V- D]>^/k 5 

f X p2 

C 

X > X 


[11 J 


where = 4. OlitlO** exp C-3464/r> 

In ■the above equations and are the free volume fraction 
of polymer and crtitical free volume fraction of polymer at which 
conversion, termination rate constant become diffusion controlled 


respectively. is calculated as follows, 

V- = -0.416 X + 0.43642 

fp 

V- = V - at X 
fx^ fp c 


[141 
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The init-ia'Lor efTiciency ratio ft is given as 

^2 


C2. ID 


where and fg are initiator efficiencies in phase 1 and phase 2 
respectively. 

2.3.2 Physical Parameters 

aD Partition coefficient for the species y between phases i and j 


fj = fiiiii] 

l[ yj]J 


C2. 2D 


e. g 


and y 


eq 

,1 


>' 12 = 


eq 


y _ _i_ 

ji 


2-1 eq 


C2. 3D 


bD Rate of mass transfer of species y from phase i to phase j 


= k 


y 

Mi J 

I 

1 1 

1 1 

w 


y 

A (— - ■ 


Mi j 

iJ 1 

''J J 


C2. 4D 


Since 


[] is known, it has been replaced by 


y . . y . 
J-J J 


y^j 1- — V — ~ above equation, can be a good estimate of 


y^l^q • When [] y^j approaches its equilibrium value, the 
driving force under equilibrium is approached. 

2.3.3 Moments and Molecular Weights 

The chain length changes during polymerization are 

characterized in terms of the moments of the growing chain length 



dist.ributi on cXj^tand of the dead polymer chain length distribution 
^k’ defined by 

- fic 

V im !/• 

C2. 5D 


^ V ^ V 

O’. CW = ) n R C13 ; Ct3 = 7 n P Ct3 

n4 " n4 " 


The number and weight average molecular weights, and M^,of 
dead polymer up to any time t are given by 


M = 
n X 


MW 


O 


m 


M = 
w X., 


MW 


m 


ca. 6D 


M 


Polydispersi ty index , PDI = 


w 




n 


C2. 7:> 


where MW is the monomer molecular weight, 
m ^ 

2.3.4 Mole Balance Equations 

Transformations of the mass balance for the various 

components involved in the polymerization using the generating 

function approach results in the following system of differential 

equations which characterize the process, 
dl ^ r ^ .1 

eft 


di . ^ r I • y ■ I • 1 

eftr " "^di ^i".I ~ "4:~' J’^Mij ^ij ^li 

J-1 '•1 J y 


i=l ,2,3 
C2.8:> 


where is the flow rate of initiator into phase i 


dM. 

1 

dtT" 


Ck . + k. D 
pi trMi ' 


M. <y-. 5 r M, M. ^ 

'^i '^i ''j J 


Mi 


i =1 ,2, . . ,5 


where F, 


Mi 


is the flow rate of monomer into phase i 


C2 . 9D 
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do- . 

af.k I. 

dt 1 di 1 


1 2 
k . . a . . 
ti Oi 

V. 

1 


3,0 


R 

r, . O'. 


-I ( ^ ^ ^ 


i=l,2,3 

C2. lOO 


do . 

;;n-^ = 2f.k..I. + k. M.C 

dt 1 di 1 trMi 1 




pi i V. 


tl 


o o 

li Oi 

V. 

1 


X X 
Oi li 


’trpi V, 


+ k 


X o 
li Oi 


trpi 


V. 

i 


i=l ,2,3 


C2. 110 


do a^. 0-. C2o. .+O-.0 

= 2f.k I. + k M.C ^ 0 + k .M. 

dt 1 di i trMi i V. pi i V, 

1 * ^ 


I. ^'Si ^Oi ,, ^Oi ^2i 

^ ^ ' « » 1C j 


ti V. 

1 


"trpi V. 


^Oi ^2i 

"trpi 


a r O'o. rV, 


dX 


Oi 


dt 




= k. ,,, M. — Si- -^k, .M. 


trMi i V. 

1 


ti i V. 
1 


w "^01 ^1. w ^ 

v.Ck. M. — 17— +k.,M^ — TT— 0 
1 trMi 1 V. tl 1 V. 

1 1 


i=l ,2,3 

C2.120 


for i=2,3 C2.130 


V V 

2 3 

where — 4~\r"0 *"*3 dead polymer 

2 3 2 3 


in phase 1. As soon as radicals gets terminated in phase 1, it 
precipitates into phase 2 and phase 3 hased on the volume 
fraction of each phase. 


*^^li _ w ^ ^ ^ 

dt trMi i V. ti V, 'trpi 


X_ , o. . 
Oi li 
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hi "oi 


O'. 


't,rpi v; "IT 


1 


^01^11 
'tl V. 


+ k 


^01^11 


11 


trpl V, 


trpl V, 


01 1 
1 J 


for i=a,3 C2. 14D 


*^^21 

di 


k, ... M. ry 

trMi i V 


2i 




°'0i °'2i 

V. 

1 


+ k 


^Oi "2i 


■trpi' 


- k 


’'21 "Oi 


t-rpi V 


Oi ^ f 

r- I 


O’. 


k. .„M. 


21 


trMl 1 V. 


+ k 


^01*^21 
tl V. 


+ k 




trpl 


’'21" 


‘t-rpl V, 


Ol 1 

1 J 


for i=2,3 C2. 15:> 


2* 3. 5 Volume change 


For Ihe monomer, polymer and skin phases ; 
dV. , . dM. . do. . dK . 


1 _ r 1 i . 1 . . liO 

dt" ~ I p ^ * p ^dtr ■*'dtr 


MW 


i=1.2,3 C2.16D 


m 


For the aqueous phase ; 



L 

^m 


MW 

n\ 




sw 


q. 5 

u w 


C2. 175 


where d and q. are volumetric flow rate of seal and injection 
sw iw '' 

water respectively. 

For vapor phase, by conservation of total volume we get, 




4 dV. 


I 


dt 


C2. 185 


2* 3 6 Conversion and Rate 

Conversion of monomer into polymer can be estimated from the 
following differential equation 
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dx 

dt 


5 dM. 


'Jo 


C2. 19:3 


Tho rate of pol ymer i zati on , r can be calculated from the 

P 

following equation , 

dx 


^ MW 
dt O m 


C2. 203 

where is the total moles of monomer charged. 

2.4 Molecular Weight Development 

The phasewise instantaneous 'number and weight average 
molecular weights at any time are calculated as follows, 

X. 


R 


ni 


Ti 


Oi 


MW 


m 


i=2,3 


C2. 213 


— ^2i 

M . = ^ MW 

wi X. . m 
li 


W, = X-, M . 
1 Oi Wi 


W, 


w, = 


i C Wg.W33 


i=2,3 

i=2,3 

i=2,3 


C2. 223 


C2. 233 


C2. 243 


Since there is no dead polymer in the phase 1, M^^ and are 

negligible. 

The overall instantaneous number average molecular weight is 


given by 


1 . O 


R 


hi 




C2. 253 


) 


The overall instantaneous weight average molecular weight M^j is 
given by 


”wl = '"2 'Jw2 * ”3 “w3= 


C2.263 


The cumulative number average molecular weight M^^is given by 
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X 


M 


_ p dx 1 .. 


where is “the terminal conversion. 

Si mi 1 ar 1 y for wei ght aver age mol ec ul ar wei ght M 


wA 


CB. 27D 


if r d:>c 1 


C2. 28^ 


2.5 Int.erphase Constitutive Relations 

Mass transfer coefficients have been calculated from 
Sherwood number estimates according to the following definition 


~MiJ 


dT. 

ij 

when Sh = 2. 0 then 


= Sh 


C2. 2o:> 


k 


y 

Mi j 


£.0^J 


or 



where d^ and are diameter of droplet and radius of particle. 

The relationship between the mass transfer coefficients k^, . and 

Ml j 

k^.. has been found out as 1271 
Mji 


,,y ^ y 

'^Mij ^Ji '^Mji 


,C2. 30!) 

These relations have been used mainly for estimating mass 
transfer coefficients. 


The mass transfer areas A. . are taken to be 




1/3 


"^12 

= C36 n 


C2. 3i:> 

^13 

= ^34 = 

n n 

P 

C2.32> 

^45 

= n R^ 

A 


C2.33D 


A^g represents the mass transfer area for the transfer between 
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monomer phase and polymer phase. represents the mass transfer- 

area for the transfer between VCM droplet surface CskinD and 
aqueous phase. Similarly the areas are for the transfer 

between monomer phase and skin phase and between aqueous and 
vapor phase r especti vel y . 


N 

i s the number of 

primary 

parti cles 

per 

uni t vol ume 

of 

initial 

monomer and this 

quanti ty 

decreases 

with 

conversion. 

As 


the polymerization temperature increases, the size of the primary 
particles at a given conversion also increases [241. At higher 
polymerization temperatures there are fewer primaries. At all 
pol ymer i zati on temper atur es the di amet er of the pr i mar i es tends 
to the same limiting value at higher conversions, approximately 
1 . 4fjm. Above a conversion of 0.03, the particles flocculate to 
produce larger units. This flocculation process can be modelled 
by 

N = N + - N C2. 345 

00 O <c 

Here t-he initial concentration is proportional to initiation 

rate and is determined mainly by temperature and 

nature/concentration of the initiator [281. The final 

concentration has been found to vary with reaction temperature 

14 3 

but it is around 2.5*10 /dm . The flocculation coefficient ^ is 
determined largely by the temperature of the process and the 
nature of the agitation. The coefficient f controls the rate of 
mass transfer of growing chai ns , monomer and initiator between 
monomer and polymer phase through the area 
2.6 Energy Balance Equations 

For the non i sothermal polymerization of vinyl chloride in a 
batch reactor an overall energy balance around the reactor yields 
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■t-h© following dynamic equa'Lions. 

The energy balance on cooling water in jacket coil is given 
below. Since coolant temperature rise is significant and flow is 


like plug flow inside the coil, an average jacket temperature T 


ja 


may be used. 


T. . + T. 
jai ' jae 


where T. . and T. are coolant inlet and exit temperatures 
respectively. 

dT . 

V,-.. P C . CT,_,- + h,.A,_CT_, - T,.> 


ja "^w w dt 


ja w jai jae 


ja ja ml ja 

C 2. 353 


Similarly for energy balance on cooling water in baffles. 


T + T 
bfi bfe 

2 

where and are coolant inlet and exit temperatures 

respectively, 

V ^w ^w dt^ = ^.fV^bfi- Tfc.^e^ " ^bfV^^m2- ^bf" 

C2. 363 

For metal wall energy balance, the temperature of the wall is 
treated as a lumped parameter. 


dT 


PmX^'K 


ml 


r . f. ^ = h. A. CT - T p + h, A, CT T. 3 

ml '^mw m dt in ja ml Ja ja ml ja 


C2. 373 


^ = h. A. ^CT - T O + 

m2 '^mw m dt in^f 




m2^ " »'bf^.f^'^m2- '^bf^ 

^ C2. 383 


where eqn C2. 373 and eqn C2. 383 are for the energy balance for 
jacket metal wall and baffle metal wall respectively. 
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For tot-al reaction mass, perfect mixing and constant holdup and 
constant physical properties are assumed in the reactor. This 
gives the reactor energy balance as 


Vo P. C 


dT 


R '"^r p dt 


= h, _A,_CT - ^ 


in ja 


+ C-AHDM, 


ml' 


m2 


dx 
O dt 


H 


C2. 395 


where H is the heat developed by agitation. 

o. 


2. 7 Branching Characteristics 

The nature and mechanism of formation of so-called anomalous 
structures in PVC have been the subject of extensive studies over 
the years. The progress made in understanding the nature of chain 
transfer to monomer and of the formation of anomalous structures 
in PVC is to a large extent due to the rapid development of new 
sensitive NMR techniques. The now accepted mechanism of chain 
transfer to the monomer has appeared in a recent paper by 
Hjertberg and Sorvik [2Q]. The reaction scheme is given in Fig 
2. 4a and Fig 2. 4b. 

According to this mechanism, the chain transfer to monomer 
starts with a head-to-head addition to the growing radical giving 
2. This new radical can take part in a propagation reaction which 
results in the formation of internal head-to-head arrangements in 
the polymer. The radical 2 is very reactive and can stabilise 
itself by rearranging to* 3. Propagation of this radical leads to 
the formation of chloromethyl side chain branching. The 
concentration of chloromethyl branches in technical PVC is about 
4-5 per 1000 monomer units. By expulsion of chlorine radical 3 
gives a polymer with a l-chloro-2-alkene end group. The chain 
transfer to monomer is completed when the chlorine atom adds a 
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CHg - CH + CHg = CH 
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Fig 2.4a Reaction Schema for chain transfer to monomer and for 
the formation of inq^ortant anomaloTis structures in PVC 
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Fig 2.4b Reaction Scheme for chain transfer to polymer by 
chlorine atom and the formation of long chain branches and 
internal unsaturation in PVC 
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monomer unit, and st-art.s a new growing chain, which will have a 
1 ,P. dichloro alkane end group C7D. Structures 5 and 7 should be 
the most frequent end groups in PVC. In ordinary PVC, the content 
of these structures is somewhat less than one per molecule. The 
content of internal head-to-head arrangements in PVC is very low 
and not higher than 0.2 per 1000 monomer units. 

Chlorine radicals can also take part in chain transfer to 
polymer. This is illustrated by the reaction scheme of Fig 2.4a. 
The first reaction in this scheme gives rise to the formation of 
on© moloc-ul© of hydrochloric acid. Attack by the chlorine atom on 
the methylene group with subsequent expulsion of chlorine offers 
an explanation of the formation of internal unsaturation in PVC 
shown by structure 12. Propagation of intermediate radicals 8 and 
9 leads to the formation of long chain branches C 0.5 per 1000 
monomer unitsD. These propagation will give long chain branches 
associated with a tertiary chlorine Cstructure 10Z> or a tertiary 
hydrogen Cstructure IID. The attack by chlorine on the methylene 
group seems to be the pref erred route. 

Long chain branches in PVC may also arise by attack of a 
macroradical on the polymer. In this case the attack is primarily 
on the chloromethyl group, giving structure 8. Other anomalous 
structures in PVC can be found by backbiting. This gives 2,4-di 
chlorobutyl branches C about 1 per 1000 monomer unitsD. The amount 
of both long and short chain branches increases as the monomer 
concentration in the polymer phase decreases. 

In spite of progress made in explaining the formation of 
short and long chain branches in PVC, the now accepted mechani sm 
of transfer to monomer and polymer are very difficult to 
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incorporate in a quantitative mathematical analysis of the 
polymerization due to the large number of unknown kinetic rate 
constants involved. In the present analysis it is assumed [141 
that transfer to polymer and i ntramolecual r_ reactions are 
responsible for the formation of short and long chain branches 
respectively in PVC. This considerably simplifies the 
quantitative analysis of the polymerization since only two lumped 
rate constants k. and k, need to be identified. Accordingly, L 
the average number of long branches per polymer chain and the 
average number of short branches per chain can be obtained from 
the solutions of below equations. 


dCX L 3 
o n 


= M Ck. /k 5 C-j— — :> ^ 
o trp p 1 - X dt 


C2. 40D 



C2.41D 


The corresponding Long Chain Branches CLCBD and Short Chain 

Branches CSCBZ) per 1000 monomer units be given by 

LCB = 1000 CL 5 C2. 42D 

n u 1 

SCB = 1000 CS X-/XO C2. 43D 

n L> 1 
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CHAPTER 3 
OPTIMIZATION 


3. 1 Objective 

The objective of the problem is to utilise the maximum 
cooling capacity of the reactor. This can be achieved by keeping 
the rate of polymerization uniform at the maximum rate possible. 
This is possible by adjusting the temperature during the 
process, according to an optimal temperature profile. 

3.2 Single Objective Optimization - 

The single objective optimization technique is based on 
Pontryagin's Continuous Minimum Principle. 

We can formulate the single objective optimization problem 


for 


Ob j ecti ve 


t 

Min I = f FCy, ui) dt 
O'* 


C3. i:> 


subject to the following constraints 

state equations : dy _ — _ y ro 

^ - y - 1 C3. 23 

control variable constraints : u . < u < u 

min max 

Let us assume there are n state variables and m control 
variables. 

The Hamiltonian is set up as 


n 


H = F + 


y A. f. 

Z i 1 


C3. 33 


i=l 


X’s are time dependent lagrange multipliers Cadjoint variables3 


The adjoint variables are governed by 

dx. _ an ■ -I o 

g^i - - i-1,2.. . . . 


,n 


C3. 45 


i=l ,2 
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Let time be floating and one of the state variable say is 

f i xed • then ter mi nal condi ti ons f or the ad j oi nt equati ons ar e 
calculated as follows, 

X^Ct^D =0.0 i=l ,2 ,Cl-kD C3. 5D 

F + Xj^fj^= O. O C3. 61) 

X^Ct^D =0.0 i=Ck+i:> n C3. 75 

The control variables are then updated for Cj+15^^ iteration as 
f ol 1 ows , 

i +1 1 

uV = u^ + c. • 3 — at any time t C3. 85 

1 1 lou. 

i = 1 , 2 m 

£ is a constant value which is computed using a sinq>lex technique 
which minimizes the objective function. 1 ^ 

3.3 Formulation for PVC Suspension Pol 3 nnnerization 

For the present optimization studies we focus on keeping the 
rate profile flat. 

Min I = fC r - r .5^ ^ dt Ca5 

qJ p d dt 


Comparing the above integrand with eqn C3. 15 

c y- ' 1 2 dx 

F = C r - r ,5 37 - 

p d dt 


Cb5 


C3. 05 


where r^ is the rate predicted and r^ is the desired rate of 
polymerization. In this case time is floating and .terminal 
conversion is fixed. Reaction mass temperature is the control 
variable. 

u = T total number of control variables m = 1 

The state variables are represented by 

y= ’®'02’®’03’^11 *°'l2'^13’*^21 '^ 22 , 

, ^ 

®23 * ^02 ' ^03 * H 2 * ^1 3*^22 ’ ^23 ’ '^l ’ ^^2 ’ ^^3 ' ^^4 * ^ 5 * ^n^O^ * ^n^ 

C3.105 
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Table 3*1 State Equations for the Optimization Problem 


The state equations are given as. 


dt'^ 

= 

dl . 

dt^ = 

f . 

J 

i=l ,a,3 

j=l .2,3 

dt'^ 

- • = 

dM. 

dt^ 

f . 

J 

i =1 , 2 * . . . ,5 

J “4 , S , • . • , < 

dt'J 

= 

d<y- . 
dt°^ ‘ 

f . 

J 

i=l ,a,3 

j=9.10,ll 

dt*^ 

=r 

^li = 
dt ^ 

f . 

J 

i=l ,2.3 

j=12.13.14 

dt'^ 

= 

do' . _ 

dT^^ 

f . 

J 

i =1 , 2 , 3 

j=lB,16.17 

dt-J 

= 

dt°‘ = 

f . 

J 

i=2,3 

j=18,19 

dt'^ 

= 


f . 

J 

i=2,3 

j=20.21 

dt'^ 

= 


f , 

J 

i=2,3 

j=22,23 

dt'’ 

= 

dV. 

dt^ 

f . 

J 

i=l ,2 5 

j=24 

dt'' 

= 

dx _ 

dt 

f . 

J 


j=29 

^J 

dt'^ 


dCX-L 3 
0 n 

= f 


j=30 


dt 

JL * 

J 



dt*^ 

= 

dS 

at^ 

f . 

J 


J=31 




4 - 4 - 














CHAPTER -4 

RESULTS AND DISCUSSIONS 

In the first section of our discussion in this chapter we 
analyse the model prediction. In the second section we have used 
our model to optimize the reactor performance using the 
formulation in Chapter 3. And in the last section we have used 
the energy balance equations reorganised in Appendix II to 
estimate the various heat transfer coefficients. 

4.1 Simulation Results of the Model 

The system of coupled ordinary differential equations 
C2. 8-2. 10, 2.40, 2. 41D have been solved numerically by the Gear 

algorithm C using NAG routine D02EBFD on HP 9000 computer. 

IL is assumed that there is no initiator available in water 
and vapor phases. Though it is true that little polymerization 
tl91 does take' place in water phase, in this present analysis it 
is assumed that there will be no polymerization in water as well 
as vapor phase also. It is also assumed that there is no transfer 
of initiator and radicals between skin and other phases. Monomer 
in water and vapor phases, though not involved in any reaction, 
is exchanged between these two phases and the other neighbouring 
phases. Moreover, the transfer of monomer to polymer phase from 
vapor phase after the critical conversion is important since it 
is accompanied by a pressure drop in the reactor. It is assumed 
that there is no dead polymer in monomer phase. As soon as 
radical gets terminated in the monomer phase, it precipitates out 
into the polymer and skin phases. 

The simulation is done under non isothermal condition. We 
used the reactor temperature profile reported by Industry for our 
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simulat/ion purpose. In t.his present work momomer charge is added 

in one slot i.e. , F.. = 0.0. 

M 

Effect of Kinetic Rate Constants 

The various rate constants used for simulation in our study 
have been reported in chapter 2. For the skin phase the rate 
constants have been taken as same in polymer phase. The initiator 
efficiency factor fg in polymer phase and f^ in skin phase are 
taken as same where as f^ in monomer phase is much higher than 
that in other phases. Although Hand. elec [101 has reported that 
kinetic behaviour of suspension PVC would change after critical 
conversion when we used his kinetic rate constants for our model 
wo did not obtain acceptable results. The polymerization rate 
decreased drastically consequently it took long time to increase 
conversion. The result of this work is shown in Fig 4.1. Xie [111 
has used hi s equati ons and f ound that Hand, el ec model 1 i ng attempts 
were premature due to the inadequate amount of kinetic data at 
these high conversion levels. 

Diff usivities and Mass Transfer Coefficients 

The necessary data related to mass transfer of species 

between the various phases C including the partition coef f icientsD 

I M M M 

are reported in Table 4.1. Diffusi vities D];..., D” , D™., D”_ are 

12 13 34 45 

estimated using standard correlations which have been referred to 
M 

in [271. is obtained from reference [181. The non - fickian 

M 

diffusion effect of D^g has been quantised and found to be one 
order less than that of fickian diffusion value. In that range we 
conducted sensitivity analysis and found it was i nsensi ti ve to 
both polymerization rate and pressure profile. The calculation of 
non - fickian diffusi vity is shown in Appendix III. 
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c o»->v'er“ siorv 

Fig 4.1 Effect of rate coiistants 
aD Sidiropolou rate constants 
b5 Hand elec rate constants 
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Table 4.1 

Necessary Data for Mass Transfer of Species Between Various 
Phases and Partition Coefficients 


Di f f usi vi ti es 

a. 

b. 


= 6. 1398wl0^‘* 

1 


^12 2 . 0»<10 


c. = 4.435W10 ® 

d. D^_ = 5. 4903»H0“® 

• 4-0 

~R _ 4 . - -1 4 

e. =10 -10 


3. Mass transfer coefficient^ 

a r 

M12 12 p 12 I 


dm^/sec 

dm^/sec 

dm^/sec 

dm^/sec 

dm^/sec 


4 fT N y 
3V, 


"]■ 


1/3 


dm/sec 


Ref 

1271 

E27] 

1273 

1271 

171 


2. Mass 

transfer areas 



a. 

*12 = *21 = tse „ M 

dm^ 

[71 

b. 

A = A = A = A 

^13 ^31 34 ^43 

dm^ 

[71 


based on droplet dia of 50/jfm and 



14 

total droplet number 4«10 



c. 

^5 = -^4 = 

dm^ 

[71 


171 


b. 


, M _ ^M r ” N V, 
M12 12 [ 3V. 


"] 


1/3 


dm/sec 


171 


R ^R r ^ N V^^ 


1/3 


dm/sec 


[71 


d. 

e. 

f . 


M 

Ml 3 
M 

M45 

y 

Mji 


“Ikt = 1-774-10-* 

l.OOOwlO”^ 

jy uy 


nj ’^Mij 


dm/sec 

dm/sec 


[271 


[271 




49 


4 . Par-tit-ion Coef f icien'ts 


a. 


= 

3. O - 5. 0 


b. 

CO 

= 

3. 04 

1271 

c . 

M 

^34 

rr 

23. 59 

C271 

d. 

M 

^46 

s: 

O. 58 

1271 

e. 

R 

^12 

= 

0, 0001 

[271 

f. 

I 

^12 

z= 

1.0 - 20.0 

[71 
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The mass transfer coefficients reported in Table 4.1 are 
obtained from the relationship Sh = 2.0» except the mass transfer 
coefficient This mass transfer coefficient is found to be 

very sensitive to pressure drop in the reactor which is discussed 
later in detail- Correlations of Ranz-Marshall » Cal der bank -Jones » 
Brian—Hales and Sh = 2.0 have been used [271 to calculate the 
mass transfer coefficients- Although the mass transfer 
coefficients calculated by different correlations are different 
the rate of pol 3 nneri zati on and pressure have been found to be 
insensitive to these coefficients in that domain. 

Partition Coefficients 

The partition coefficients have been reported in Table 4.1. 

M I 

Among these partition coefficients, have been found 

to be sensitive to polymerization rate, this is discussed later 
i n detai 1 . 

As noted in section 2.2.1, at conversions less than 0.001, 
the polymerization takes place only in monomer rich phase. So, up 
to 0.001 conversion only monomer phase equations are simulated 
and after that the entire model is included. The skin phase is 
assumed to be formed at 0.001 conversion with the initial 
thickness of O, l/um i . e. , a mono particle layer of polymer tlOl. 

The physical parameters used for simulation has been 
reported in Table 4.2. 

Sensitivity Analysis 

The effect of various model parameters on rate of 
polymerization and pressure have been studied to find out which 
parameters are more sensitive. 
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Table 4,2 

Physical Parameters Used in Simulation 


Parameters 

it 

Val ues 

Uni ts 

Ref 

1 . 

A, -/A . 
bf ja 

0. 28940 


[31 1 

2. 

C 

P 

3. 32 

J/gm. K 

[271 

3. 


o 

CD 

J-/gm. K 

[27J 

4. 

S 

4. 18 

J/gm. K 

[27J 

5. 


0. 0171 

-3 


[ 31 1 

6. 

V yy 
m2 R 

2-0»H0 


[311 

7. 

V /V 
bf R 

0. 01427 


131] 

8. 

V. /V- 
ja R 

0. 0312 


[31 1 

Q. 

^m 

840.2 

gm/dm^ 

[91 

10. 


1370. 3 

gm/dm^ 

[Q] 

11. 

^mw 

7850. 0 

gmx'dm^ 

[271 

12. 


085. 44 

gmx'dm^ 

[01 

13. 


62. 45 

gm/gmol e 

[01 

14. 

MWj 

346. 45 

gm/gmol e 

[311 

15 

H/Dr 

1 . 2934 


[311 


C£ 'i 'l TK ■ L L ! BRA R V 

! ’ " . KANPUR 


^jvo. A. J1Q2G0 
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The object-j ve function for this sensitivity analysis has 
been f ormul ated as foil ows > 


I = f Cr - r ^ dt 
1 O'* P ^ 


^3 = CPp - Pd5^'° dt 


where and are "the integral errors in polymerization rate 
and pressure respectively; r and r , are the polyinerization rate 

P Q 

predicted by simulation and reference rate profile. Similarly p^ 
and p^ are the pressure predicted by simulation and reference 
pressure profile. 

For each parameter lower and upper bound values have been 
found out by simulation experience. With these lower and upper 
bound values for each parameter > the polymerization rate and 
pressure have been simulated and the integral error and Ig are 
calculated. 

The various model parameters like kinetic parameter fty 

diffusivities mass transfer coefficients ’'M45M 

and partition coefficients ^34* ^45’ ^12 been 

M 

taken for this analysis. We found among these parameters, 

I R 

y^g, ft are more sensitive to polymerization rate Cin the 

same order as indicated!) but insensitive to pressure profile. On 

M 

the other hand has been found to be sensitive to pressure 

profile but insensitive to polymerization rate. All other 
parameters are neither sensitive to polymerization rate nor to 
pressure profile. 


Effect of y. 


12 


The monomer pjartition coefficient between phase 1 and 
M 

phase 2 , Is the most sensitive parameter to rate of 
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polymerization. A slight increase in the monomer concentration of 

polymer phase is enough to accelerate the reaction .The monomer 

concentration in the polymer phase is very important in 

determining the rate curve. The rate decreases drastically with 

M 

increase in the value of The effect of this on rate is shown 

in Fig 4.2. Since monomer phase disappears after critical 
conversion C x^= 75% > the monomer in the polymer phase starts 
depleting thereby producing a steep fall in rate. 

Effect of 7 ^^ 2 ; 

As increases the rate of polymerization decreases. 

If high then less amount of initiator is located in the 

polymer phase » therefore rate decreases where as for low value 

the initiator concentration in polymer phase is high so rate 

I 

increases. Although the rate profile does vary with different ^ 
values, i't is not, very significant,. The reason for this is the 

initiator location is irrelevant for high radical diffusivity 

j? —0 2 I 

C = 2. 0»tl0 dm /secD . The effect of is shown in Fig 4.3. 

Effect of Dj g 

The rate of radical mass transfer is crucial in determining 

R 

the shape of the reaction rate profile. For a high value of 
in the range 10 ~ 10 radicals initiated in the monomer phase 

1 diffuse rapidly into the polymer phase 2 and so essentially 

all growth takes place in phase 2. On the other hand for low 

R —10 -14 2 

radical diffusivity C in the range lO - 10 dm /sec 0 

essentially there is no radical transfer, thus the rate peak is 

R 

lowered. The effect of on polymerization rate at different 

' R 

values is shown in Fig 4.4. We found that has negl 1 gi bl e 

■—5' "'—S''"'- 

effect on rate of polymerization in the range 10 - 10 . But in 
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—8 “ 1 4 

the range 10 - 10 the rate profile decreases with decrease 

in 

Effect of Initiator Efficiency Ratio ,/? 

The increase in initiator efficiency ratio has been found 
to increase the polymerization rate. As (3 increases more amount 
of initiator in the polymer phase decomposes to ^yield more 
radicals which triggers the polymerization reaction. This effect 
on rate is shown in Fig 4.5. 


Effect of 

The monomer interphase mass transfer coefficient between 

aqueous and vapor phase plays a vital role in determining the 

M 

pressure drop in the reactor. If is very low then transfer 

of monomer into vapor phase is negligible thereby causing 

pressure to drop faster. On the other hand if it is high more 

amount goes into vapor phase reducing the pressure drop. Fig 4.6 

M 

shows the effect of on pressure profile. 

Parameter Estimation 

The objective function for parameter estimation has been 
formulated as follows. 


Min 


''f 

= ECr 
i=l ^ 




2.0 


where r^ and r^ are the average rate of polymerization predicted 
by simulation and average rate of polymerization given by 
industry at every hour Cup to t^^ hour>. 

With initial estimates C Table 4. 3;> for the rate sensitive 
parameters the rate of polymerization is simulated and averaged 
over every hour. With this the error is calculated. The 








corw/ersfon 


Fig 4.4 Effect of radical diffusivity 

a!) D?- = 1.05H0~® : b> = 1.0*1 o“® ; c3 = 1.0*10 








con V'er Sion 


Fig -4. 6 Effect. of monomer Interphase inas 5 s 

M 

coefficient, 

M4o 

a:> = l.Oitlo"® ; W = 1.0X10'^^ 


transfer 
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Table 4.3 Model parameter values used before and after 
parameter tuning 


Parameters 

Befor e Tuni ng 

Af ter Tuni ng 

Ref 


6. 1308»fl0^‘^ 

6. 1398x10^'^ 

1271 


—1 "5 

3. 0*10 

—1 "5 

3.0x10 

1271 


4.435*«10“® 

4. 435x1 0“® 

1271 


5. 4903x1 0~® 

5. 4903x10"® 

£271 


4. 435x10"® 

4.435x10"® 

£271 

-fa 

1. Oxio"*^ 

3.0410"® 


M 

^13 

3. 1 

3.31 


M 

^13 

3. 04 

3. 04 

£271 

M 

^34 

33. 59 

33. 59 

£271 

M 

^45 

0. 58 

0. 58 

£271 

I 

^13 

3. 448 

3. 7831 


R 

^13 

0.0001 

0. OOOl 

£271 

ft 

0.5 

0.65 


M46 

l.Oxlo"'^ 

1.707440"® 






usi ng si mpl ex 


parame-ters are then tuned, to minimize 
opti mi zati on techni que , 

Similarly pressure sensitive parameters are tuned to match 
simulated pressure with industrial pressure data. 

The parameter estimation has been done for seven industrial 
runs. In all the runs the parameter values turn out to be nearly 
same, allowing us to take an average value for each parameter. 
The parameter values for seven runs is given in Table 4.4 

With the estimated parameters the results are simulated. 
These results like rate of polymerization, pressure drop in the 
reactor, and conversion history matches well with industrial 
data. The comparison of industrial and our simulated conversion 
history, polymerization rate and pressure profile are shown in 
Fig 4. Fig 4. 8cL^and Fig 4. QE* For second run the respective 
results are shown in Fig 4.7b, Fig 4.8b and Fig 4.9b. 

4*2 Optimization Results 

The system of coupled ODE’s C Table 3. ID have been solved 
using NAG routine E>02EBF. The derivative in equation C3. 16D 

has been obtained numerically using NAG routine EX)4AAF. Again the 
system of coupled ODE’s equation C 3.1 2D has been solved using 
D02EBF routine. The weighting factor e has been calculated using 
simplex technique E04CCF in NAG. 

First with the initially assumed temperature profile to be 
flat at 330K the equations in Table 3.1 have been solved to get 
the state variable values. These equations are solved till 
conversion reaches 80%. The reason for choosing terminal 
conversion as 80% is above this conversion it is not fx3®si bl e to 
keep the rate profile flat as there is virtually no monomer in 
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Table 4.4 

Parameter 

Estimation For 

Seven Runs 


Runs 

M 

I 

^la 

*^13 

ft 

M45M 

1. 

3. 1065 

4. 038 

a. 08a»tio“® 

0.65 

1. 70»tl0“® 

a. 

3. ai60 

3. 4073 

1. 401*<1O“® 

0. 64 

1.70*1 o“® 

3. 

3.1000 

3. 804 

5. 600»tlo“® 

O. 643 

1.70*10~® 

4. 

3. aoso 

3. 6oa 

a. ooo»fio”® 

0. 65 

1 . 70*1 0“® 

5. 

3. 1016 

4. 378 

8. 000»H0~^ 

O. 646 

1.70*10“® 

6. 

3. aaol 

4. 038 

a. 68Q»10”® 

0.640 

1.70*10“® 

7. 

3. ai43 

3. 4167 

5. 000»tio“^ 

O. 653 

1.-70*10“® 




•time <miri) 


Fig A. 7a Conqjarison of conversion history between 
experimental data and model prediction for run 1 
n, e>qDeri mental data ; 


model 
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•time <mlr») 

^i-9 4-. 7b Conq>arison of conversion hist-ory beiween 

experimental data and model prediction for run 2 
o» experin^ntal data : - model 



•time <min) 


Fig 4# 8a Comparison of polymerizalion rate between 
experi mental data and model prediction for run 1 
o, experimental data ; model » 











■time Cmln) 

Fig 4.9b Comparison of pressure profile between 
experimental data and model prediction for run 2 
a, experimental data ; model 
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the polymer phase. After solving the state equations forward, the 
Hamiltonian is set up as given in Chapter 3 and the adjoint 
equations C3. i2D are integrated backwards. Then, the temperature 
profile is updated using equation C3. 16>. To attain the 
convergence quickly optimal value of e is calculated at every 
i ter ati on by si mpl ex techni que. 

Optimization Result vs Industrial Result 

The optimized rate profile has been found to be fairly flat 
which is shown in Fig 4.10. In this case the heat release is 
uniform and at the maximum cooling capacity. By operating the 
reactor along optimal temperature profile a 10% saving in 
reaction time can be achieved. The comparison of two conversion 
histories are shown in Fig 4.11. The corresponding reactor 
operating temperatures are shown in Fig 4.12. 

Molecular Weight Averages 

The cumulative molecular weight averages for the optimal 
temper atur e pol icy condi ti on are 1 i t tl e 1 ess than that of 
industrial run. This is shown in Fig 4.13. The corresponding 
polydispersi ty index CPDI3 is shown in Fig 4.14. It has been 
found that for industrial run PDI is around 2. 05 where as for the 
optimal temperature condition it is around 2.1, 

Short and Long Chain Branching Characteristics 

The formation of long chain branches CLCB3 per thousand 
monomer units both for the industrial and optimal temperature 
policy run have been found to be very negligible. On the other 
hand, the formation of short chain branches CSCB> is appreciable. 
For the optimal temperature policy run it is slightly higher than 
for industrial case. This is shown in Fig 4.15. 









Cumulat-iva molecular weight, averages 
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Fig 4*13 Comparison ot cumulative molecular weight averages 
a!) industrial molecular weight averages 

bZ) optimal temperature policy molecular weight averages 




con\/ersion 


4^*1 ■4 Con^sarlson of polydisporsity index 
aD industrial PDI 

b;> optimal temperature policy PDI 



• of dof ocbi V© uni "ts/l 
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4*3 Heat Transfer Coefficients Results 

To estimate the various heat transfer coefficients, the 

energy balance equations in Chapter 2 have been reorganised and 

given in Appendix II. In this system of equation the unknowns are 

jacket side heat transfer coefficient h. » reactor side heat 

ja 

transfer coefficient h. , and baffle side heat transfer 

xn 

coefficient h, The known variables are T. , T, - and T. 

bf j a bf 

Closed Loop Real Time Optimization CCLRTOZ> technique 

suggests to solve the five simultaneous equations coupled with an 

optimization technique to estimate the unknown parameters 1301. 

First initial values for these parameters are assumed and then 

the residuals of the simultaneous equations are driven to zero. 

The calculated variables like T, T. and T,^ are compared with 

the experimental data. If the values does not match, the 

parameter values are updated by using an optimization technique. 

Since the coolant flow rate, inlet and outlet temperature of 

the coolant in the jacket is constant over the time we have 

2 

fixed h. at 40.0 W/dm . K. This value has been obtained by trial 
and error method. 

Since we have fixed h. we are left with two unknown 

parameters h^^^ and h^.^. The initial values for these parameters 

2 

have been assumed as 40. 0 and 20. 0 W/dm .K respectively. With 

these values the equations have been sol ved simultaneously using 

NAG routine C05NBF. The parameters are then updated using NAG 

routine E04CCf , Simplex optimization technique. 

The results are shown in Table 4. 5. The reactor inside heat 

transfer Goefficient h. thus estimated found to vary between 9.5 
2 ' 

- 10. 8 W>dm . K whose value is close to the literature value IT . 5 



W/^dm .K 1161 . h. remains uniform ‘through out the reaction time 

in 

because continuously water is added into the reactor which keeps 
the mass slurry. The agitator power consumed during reaction 
remains almost uniform as per the industrial data. This shows 
that the viscosity of the medium does not rise significantly as 
polymerization proceeds, has been found to vary with time 

because the coolant flow rate varies with time. 




Table 4.5 Heat transfer coefficients result 

IB 

TIME 

T. CO 
ja 

^bf 

T CO 


40 

33. 34000 

37. 69000 

58. 36990 


45 

33. 38000 

44. 24000 

58. 36000 


50 

33. 42000 

47. 50000 

58. 36990 


55 

33. 44000 

44. 64000 

58. 36000 


60 

33. 46000 

41 . 02000 

58. 36000 


65 

33. 50000 

38. 15000 

58. 38000 


70 

33.51000 

37. 02000 

58. 36000 


75 

33. 48000 

34. 72000 

58. 20900 


80 

33. 47000 

32. 69000 

58. 29990 


85 

33. 45000 

32. 55000 

58. 24000 


00 

33. 42000 

31 . 64000 

58. 18000 


05 

33. 38000 

30. 46000 

58.13000 


100 

33. 35000 

30. 62000 

58. 08000 


105 

33. 32000 

20. 60000 

58. 02000 


110 

33. 28000 

30. 32000 

57. 97000 


115 

33. 24000 

30. 74000 

57. 93000 


120 

33.21000 

30. 87000 

57. 84990 


125 

33. 19000 

30. 81000 

57. 79990 


130 

33.16000 

31 . 58000 

57. 75000 


135 

33.13000 

31 . 48000 

57. 70000 


140 

33.11000 

31 . 20000 

57. 63990 


145 

33.10000 

31.13000 

57. 59000 


150 

33. 07000 

31 . 05000 

57. 52000 


155 

33. 06000 

31 . 35000 

57. 47990 


160 

33. 07000 

31 . 06000 

57.41990 


165 

33. 08000 

32. 45000 

57. 36000 


170 

33. 09000 

32. 24000 

57. 29990 


175 

33. 09000 

30.46000 

57. 22990 


180 

33. 07000 

33. 06000 

57.18000 


185 

33. 06000 

31 . 83000 

57.13990 


10O 

33. 04000 

32. 44000 

57. 09000 


105 

33. OlOOO 

32.38000 

57. 02990 


200 

32.98000 

32. 04000 

56. 97000 


205 

32. 95000 

32.11000 

56.91990 


210 

32. 92000 

32. 28000 

56. 86000 


215 

32. 89000 

32.17000 

56. 7^00 


220 

32. 87000 

32. 49000 

56.74000 


225 

32. 86000 

31 . 82000 

56. 70000 


230 

32.84000 

31 . 70000 

56. 65000 


235 

32.81000 

31 . 97000 

56.59^0 


240 

32. 78000 

31. 15000 

56.54000 


245 

32. 76000 

30. 77000 

56.47000 


250 

32.73000 

31 . 28000 

56. 43000 


255 

32.72000 

30.96000 

66.40000 


260 

32. 70000 

30.68000 

56.38990 


265 

32.67000 

2Q . 87000 

56.38000 


270 

32.63000 

28.73000 

66.36M>0 


275 

^.60000 

28.01000 

56.36990 


280 

32. 57000 

27. 38000 

56.36000 


285 

32.54000 

27. 22000 

56.34990 


290 

32.51000 

28.00000 

56.34000 


295 

32.49000 

29.87000 

56.31900 




TIME 


h 


bf 


h. 


in 


40 

4. lOOOO 

9. 53291 

45 

10.10000 

9. 86836 

50 

9. 50000 

9. 94772 

55 

16.20000 

10.04072 

60 

14. 00000 

10. 02189 

65 

21 . 00000 

9. 95099 

70 

15. OOOOO 

10.02066 

75 

27. OOOOO 

10.11341 

80 

26. OOOOO 

10.27955 

85 

26. OOOOO 

10.17785 

so 

27. OOOOO 

10.13273 

95 

28. OOOOO 

10.07891 

100 

33. OOOOO 

10. 05702 

105 

38. OOOOO 

10.03902 

110 

41 . OOOOO 

9. 98032 

115 

37. OOOOO 

9. 89230 

120 

39. OOOOO 

9. 92289 

125 

42. OOOOO 

9. 88665 

130 

40. OOOOO 

9.83645 

135 

41 . OOOOO 

9.81968 

140 

41 . OOOOO 

9. 77907 

145 

42. OOOOO 

9. 82238 

150 

42.00000 

9. 78052 

155 

42. OOOOO 

9.77681 

160 

41.00000 

9.81816 

165 

41 . OOOOO 

9.86909 

170 

42.00000 

9. 95907 

175 

41 . OOOOO 

9. 97349 

180 

40. OOOOO 

9. 90194 

185 

42. OOOOO 

9. 90917 

100 

42. OOOOO 

9. 89632 

195 

42.00000 

9.86149 

200 

42.00000 

9. 80785 

205 

40. OOOOO 

9. 76768 

210 

42. OOOOO 

9.75300 

215 

42. OOOOO 

9. 71542 

220 

42. OOOOO 

9. 67740 

225 

42. OOOOO 

9. 69553 

230 

43. OOOOO 

9. 67286 

235 

43, OOOOO 

9. 62655 

240 

43.00000 

9. 59385 

245 

43. OOOOO 

9.59194 

250 

43.00000 

9. 51848 

255 

43. OOOOO 

9. 51786 

260 

43. OOOOO 

9. 47599 

265 

43.00000 

9. 41290 

270 

43. OOOOO 

9.33206 

275 

43. OOOOO 

9.27774 

280 

43.00000 

9. 21239 

285 

43.00000 

9.13349 

290 

43.00000 

9.C©105 

295 

43. OOOOO 

8.98963 


TIME 




40 

60796. 11588 

1469. 56588 

45 

60768. 75417 

8618.73194 

50 

60741 . 66667 

8870. 69861 

55 

60850. 83750 

4016. 58375 

60 

60686. 67083 

3737. 50048 

65 

60876. 58778 

5518. 73611 

70 

60331 . 85000 

4355. 88889 

75 

60638. 88888 

7507. 86153 

80 

61184. 44861 

7685.51431 

85 

60850. 83750 

7693. 91708 

90 

60883. 75000 

8071 . 50048 

95 

60796. 11588 

8578. 83611 

100 

60883. 75000 

9510.78888 

105 

60883. 75000 

10898. 43097 

110 

60883. 75000 

11308.87500 

115 

60768. 75417 

10890. 51389 

180 

60850. 83750 

10876. 04167 

185 

60883. 75000 

11401.37500 

130 

60796. 11588 

10788. 98653 

135 

60883. 75000 

11831.70875 

140 

60768.75417 

10947. 15319 

145 

60850. 83750 

11373.98653 

150 

60883. 75000 

11889.19444 

155 

60850. 83750 

11068. 09788 

160 

60796. 11588 

10843.18097 

165 

60883. 75000 

10681.55597 

170 

60905. 83333 

11081. 85000 

175 

60850. 83750 

11880.79167 

180 

60796.11588 

10148. 83611 

185 

60768.75417 

11130. 50000 

190 

60796.11588 

10919. 79808 

195 

60883. 75000 

10758.36153 

800 

60796.11588 

10698. 16708 

805 

60796.11588 

10375.33333 

810 

60883. 75000 

10977. 85048 

815 

60883. 75000 

11086. 50048 

880 

60796. 11528 

10668. 09788 

885 

60883. 75000 

11144. 15319 

830 

60883. 75000 

11394. 66667 

835 

60823. 75000 

11861.80597 

840 

60796. 11588 

13130. 94486 

845 

60883. 75000 

11710.58819 

850 

60883.75000 

11349. 36153 

855 

60883. 75000 

11650.33375 

860 

60883. 75000 

18017. OOOOO 

865 

60796. 11588 

18848. 75048 

870 

60796.11588 

13874.79808 

875 

60850. 83750 

14104. 6^48 

880 

60850.83750 

14385.80875 

885 

60883.75000 

14844.19444 

890 

60850.83750 

18434. 65319 

895 

60850.83750 

10585.81944 
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CHAPTER 5 

CONCLUSIONS AND SUGGESTIONS 

B. 1 Conclusion 

A mulii -phase mathematical model has been presented for free 
radical suspension polymerization of VCM in a batch reactor. This 
model takes care of precipitation of monomer from monomer phase 
into other phases and predicts pressure drop in the reactor. We 
have assumed here that VCM vapor obeys ideal gas law. Apart from 
pressure prediction this model also predicts polymerization rate, 
cumulative molecular weight averages and long and short chain 
branches characteristics. 


There are several parameters in the model. The sensitivity 
of these model parameters to polymerization rate and pressure 
profile has been studied. For this purpose, by simulation 
experience we fixed the upper and lower bound values for each 
parameter and sensitivity analysis has been investigated in that 


domain. Based on our sensitivity analysis we found that four 


parameters viz. , monomer partition coefficient initiator 

I R 

partition coefficient 7 ^^^* radical diffusivity And initiator 

efficiency ratio (3 were very sensitive to polymerization rate but 


insensitive to pressure profile. On the other hand, we found 

M 

monomer interphase mass transfer coefficient k^.-? was sensitive 

M4-0 

to pressure profile but insensitive to polymerization rate/^^ ^ ^ 


Among these parameters 7 ^^^ is the most sensitive parameter 


which determines the rate profile. We found by increasing the 
monomer concentration in the polymer phase the rate could be 
enhanced . Thi s suggests a semi —batch mode of oper ati on wi th 
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respect to monomer after critical conversion to keep the monomer 
concentration in the pol 3 rmer phase as high enough to prevent a 
fall in pressure. This would also minimize the formation of 
defective structures thus improving the thermal stability of PVC. 
We could not investigate this semi -batch mode operation as we ran 
out of time. 

Having found out the sensitive parameters, we estimated 
these val ue usi ng i ndus tr i al data . The Si mpl ex opt i mi zati on 
method has been used to estimate the value of these parameters . 
We estimated the parameters for seven such industrial run and an 
average value has been taken for each parameter. 

To utilise the maximum cooling capacity of the reactor, we 
found an optimal temperature policy using Pontryagin^s Continuous 
Minimum Principle. When the reactor is operated along this 
optimal temperature policy a 10^ saving in time is realised. 
Although there is a increase in productivity this has it^s own 
drawbacks. By operating along optimal temperature policy the 
cumulative molecular weight averages decreases slightly which in 
turn may affect the grade of PVC resin. If this fall is not 
acceptable for the end use of PVC it can be incorporated as a 
constraint- 

All the three heat transfer coefficients Ch . , h,^, h. D in 

j a Di 1 n 

the energy balance have been estimated using industrial data. The 

estimated heat transfer coefficient h. lies close to the 

in 

reported value in the literature. 
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5» 2 Sugges'tions for future work 

iD To use multiple initiator systems for achieving flatter heat 
release. To utilise the maximum cooling capacity of the 
reactor it is necessary to keep the rate profile flat. The 
dominating approach to this problem is to use a combination 
of rapid and slow initiators. 

iiO Agitation is very important in the suspension PVC process 
in that it determines the monomer droplet size. It governs 
the stability of the suspension during polymerization 
and particle size of the product formed. It also plays an 
important role in determining product properties such as 

porosity and bulk density. Also, it has a significant effect 
on the reactor inside heat transfer coefficient. In general 
more agitation will increase inside heat transfer 

coefficient and, hence, the overall heat transfer 

coefficient. Therefore modelling related to agitation is of 
prime importance. 

iii> To incorporate the quantitative relationship between product 
characteristics Clike porosity, particle size distribution, 
crystallinity, intrinsic viscosity^, hydrodynamic and 

operating conditions of the batch reactor. 

iv> To operate the reactor in a semi -batch mode to maintain 

monomer concentration as high enough to prevent a drop 
in the ojjerating pressure. This drop occurs nearly at 
critical conversion. Till critical conversion batch mode is 
followed and after switched to semi -batch mode to reduce the 
formation of defective structures, prevent thermal 
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instability and to prevent fall in the inolecular weight of 
PVG. It is suggested to add monomer in the vapor phase 
continuously in the semi —batch mode so that the polymer 
particle sucks the monomer from vapor phase. 
v3 To use energy balance equation in the model to predict the 
reaction temperature. This is essential to control the 
reactor temperature along the optimal temperature profile. 
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APPENDIX I 


Calculation for initiator distribution among three phases on 
initiator addition 

Let the total flow rate of Initiator be F j . This amount gets 
distributed among three phases C monomer, polymer and skin 
phases 3 based on initiator partition coefficients between 
different phases, is taken as same as 

F=F+F+F CIID 

I II *^12 13 vx.ij 

where Fj * Fjg* are the amount of initiator flowing into 

phase 1. phase 2 and phase 3 respectively. Let At be the time 
interval for integrating the state equations C2. 8-2. 19,2. 40,2. 413 
We know , .. 




+ F. 


II 


+ F. 


12 


At 3/V^ 

At37^ 


CII.23 


vdiere I^ and I^ are the initiator amount already present in 
phase 1 and phase 2 before adding up further. 

Also, 


^13 ^12 


C I. + F_.At 3/V, 
1 11 1 


^3 * '' 13 **- 


CII.33 


Knowing Fj and we can solve the above three equations 

simultaneously to yield F^^, F^g and Fjg. I^ , Ig, Ig, Vg and 


Vg are obtained from while simulating the model. 
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APPENDIX II 


Formulations for Heat Transfer Coefficients Estimation 


AT . 

R, = m. CCT. . - T. 5 + h. A. CT. - T. :> - V. pC 

1 ja w jai jae ja ja ml ja ja'^w w At 


AT, 


■'a = " V*bf’^^»2 - ■^bi-" - ''bf''w‘=v,-^ 


bf 


AT 


R.^=h. A. CT-T„D - h.A. CT„ - T.D - p CV 
3 an ja ml ja ja ml ja m ml At 


ml 


AT 


^4 = ^n^,f<'‘^ - “^ma" - »'bf^>f^T^ " T^f^ - ^mw^m^ma-S*: 


ma 


dx 


R_ = C-AH3M- ^ - h. A. CT - T ,3 - h. A. -CT - T 
5 O dt an ja ml an of . ma 


- V p C ^ + H 
R r p At a 

where H is the heat developed by the agitation. It is estimated 

Si 

[(61 that 90^ of agitator power is converted into heat. 

This open equation based model has been formulated to suit 

CLRTO C Closed Loop Real Time C^timizationD techni que [ 301 . R^ , 

R^, and R„ are the residuals of the simultaneous 

3 3 4 5 

equations. 
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APPENDIX III 

Quanlisa-tion of non - flckian diffusion of VCM in PVC £181 







The above graph shows the fickian and non - fickian 
diffucion process of VCM in PVC. 

The fickian diffusion of VCM in PVC can be calculated as 
follows » 


Y ^ ^ a a. _a- cm . id 

= 1 - C6/fr D } Cl/n D e>q> ^ ^ 

«e n^ 


where M^ is the weight of'‘ penetrant entering or leaving the 

sphere in time t after an instantaneous change in surface 

concentration, M is the total weight change after infinite time 

at the new surface concentration, n is the series of integers, 
M 

pr. is the diffusion coefficient, and d is the diameter of 

la 

the sphere. 

For calculating fickian diffusion consider a F>oint on 

fickian diffusion curve. At t^' ® = ■4. 0min^'^, M^XM^ = 0.00. 

—5 

The diameter of particle is d = 4.47x10 dm. Substituting these 
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values in eqn CIII.IZ) we get = 2. 0*^10 dm^/sec. 

For quantising the non - fickian diffusion process, let us 
consider a point on the non - fickian curve in the graph shown 
here. 

At t®'® = 4 min® ®, M. =3.0 mg/gm , M = 3. 4 mg/gm 

t/ oc 

Substituting these values in eqn Cl II. 13 
we get = 8. 66»*10 dm^/sec. 





APPENDIX IV 


" I ' 

StatiiiVy of Conq^utat-ional Technique 

The stability of numerical technique used in opti jni zati on 


problem has been investigated by reducing the step length of time 
At from 200 sec to 100 and 50 sec. The effect of this reduction 
on optimal temperatue profile is shown below. From this graph one 
can conclude that the present step length At = 200 sec is a good 
increment in time. 



Fig IV. 1 Effect of step length At 

a> At = 200 b:> At = 100 c5 At = 50 sec. 



